The transcriptional cofactor Tle3 reciprocally controls effector and central memory CD8+ T cell fates.

Antigen-experienced CD8+ T cells form effector and central memory T cells (TEM and TCM cells, respectively); however, the mechanism(s) controlling their lineage plasticity remains incompletely understood. Here we show that the transcription cofactor Tle3 critically regulates TEM and TCM cell fates and lineage stability through dynamic redistribution in antigen-responding CD8+ T cell genome. Genetic ablation of Tle3 promoted CD8+ TCM cell formation at the expense of CD8+ TEM cells. Lineage tracing showed that Tle3-deficient CD8+ TEM cells underwent accelerated conversion into CD8+ TCM cells while retaining robust recall capacity. Tle3 acted as a coactivator for Tbet to increase chromatin opening at CD8+ TEM cell-characteristic sites and to activate CD8+ TEM cell signature gene transcription, while engaging Runx3 and Tcf1 to limit CD8+ TCM cell-characteristic molecular features. Thus, Tle3 integrated functions of multiple transcription factors to guard lineage fidelity of CD8+ TEM cells, and manipulation of Tle3 activity could favor CD8+ TCM cell production.

Antigen exposure reshapes chromatin architecture in central memory CD8+ T cells and imprints enhanced recall capacity.

CD62L+ central memory CD8+ T (TCM) cells provide enhanced protection than naive cells; however, the underlying mechanism, especially the contribution of higher-order genomic organization, remains unclear. Systematic Hi-C analyses reveal that antigen-experienced CD8+ T cells undergo extensive rewiring of chromatin interactions (ChrInt), with TCM cells harboring specific interaction hubs compared with naive CD8+ T cells, as observed at cytotoxic effector genes such as Ifng and Tbx21. TCM cells also acquire de novo CTCF (CCCTC-binding factor) binding sites, which are not only strongly associated with TCM-specific hubs but also linked to increased activities of local gene promoters and enhancers. Specific ablation of CTCF in TCM cells impairs rapid induction of genes in cytotoxic program, energy supplies, transcription, and translation by recall stimulation. Therefore, acquisition of CTCF binding and ChrInt hubs by TCM cells serves as a chromatin architectural basis for their transcriptomic dynamics in primary response and for imprinting the code of "recall readiness" against secondary challenge.

Low RNA stability signifies increased post-transcriptional regulation of cell identity genes.

Cell identity genes are distinct from other genes with respect to the epigenetic mechanisms to activate their transcription, e.g. by super-enhancers and broad H3K4me3 domains. However, it remains unclear whether their post-transcriptional regulation is also unique. We performed a systematic analysis of transcriptome-wide RNA stability in nine cell types and found that unstable transcripts were enriched in cell identity-related pathways while stable transcripts were enriched in housekeeping pathways. Joint analyses of RNA stability and chromatin state revealed significant enrichment of super-enhancers and broad H3K4me3 domains at the gene loci of unstable transcripts. Intriguingly, the RNA m6A methyltransferase, METTL3, preferentially binds to chromatin at super-enhancers, broad H3K4me3 domains and their associated genes. METTL3 binding intensity is positively correlated with RNA m6A methylation and negatively correlated with RNA stability of cell identity genes, probably due to co-transcriptional m6A modifications promoting RNA decay. Nanopore direct RNA-sequencing showed that METTL3 knockdown has a stronger effect on RNA m6A and mRNA stability for cell identity genes. Our data suggest a run-and-brake model, where cell identity genes undergo both frequent transcription and fast RNA decay to achieve precise regulation of RNA expression.

MLL3/MLL4 methyltransferase activities control early embryonic development and embryonic stem cell differentiation in a lineage-selective manner.

H3K4me1 methyltransferases MLL3 (KMT2C) and MLL4 (KMT2D) are critical for enhancer activation, cell differentiation and development. However, roles of MLL3/4 enzymatic activities and MLL3/4-mediated enhancer H3K4me1 in these processes remain unclear. Here we report that constitutive elimination of both MLL3 and MLL4 enzymatic activities prevents initiation of gastrulation and leads to early embryonic lethality in mice. However, selective elimination of MLL3/4 enzymatic activities in embryonic, but not extraembryonic, lineages leaves gastrulation largely intact. Consistent with this, embryonic stem cells (ESCs) lacking MLL3/4 enzymatic activities can differentiate toward the three embryonic germ layers but show aberrant differentiation to extraembryonic endoderm (ExEn) and trophectoderm. The failure in ExEn differentiation can be attributed to markedly reduced enhancer-binding of the lineage-determining transcription factor GATA6. Furthermore, we show that MLL3/4-catalyzed H3K4me1 is largely dispensable for enhancer activation during ESC differentiation. Together, our findings suggest a lineage-selective, but enhancer activation-independent, role of MLL3/4 methyltransferase activities in early embryonic development and ESC differentiation.

Abnormalities in intron retention characterize patients with systemic lupus erythematosus.

Regulation of intron retention (IR), a form of alternative splicing, is a newly recognized checkpoint in gene expression. Since there are numerous abnormalities in gene expression in the prototypic autoimmune disease systemic lupus erythematosus (SLE), we sought to determine whether IR was intact in patients with this disease. We, therefore, studied global gene expression and IR patterns of lymphocytes in SLE patients. We analyzed RNA-seq data from peripheral blood T cell samples from 14 patients suffering from systemic lupus erythematosus (SLE) and 4 healthy controls and a second, independent data set of RNA-seq data from B cells from16 SLE patients and 4 healthy controls. We identified intron retention levels from 26,372 well annotated genes as well as differential gene expression and tested for differences between cases and controls using unbiased hierarchical clustering and principal component analysis. We followed with gene-disease enrichment analysis and gene-ontology enrichment analysis. Finally, we then tested for significant differences in intron retention between cases and controls both globally and with respect to specific genes. Overall decreased IR was found in T cells from one cohort and B cells from another cohort of patients with SLE and was associated with increased expression of numerous genes, including those encoding spliceosome components. Different introns within the same gene displayed both up- and down-regulated retention profiles indicating a complex regulatory mechanism. These results indicate that decreased IR in immune cells is characteristic of patients with active SLE and may contribute to the abnormal expression of specific genes in this autoimmune disease.

CTCF mediates CD8+ effector differentiation through dynamic redistribution and genomic reorganization.

Differentiation of effector CD8+ T cells is instructed by stably and dynamically expressed transcription regulators. Here we show that naive-to-effector differentiation was accompanied by dynamic CTCF redistribution and extensive chromatin architectural changes. Upon CD8+ T cell activation, CTCF acquired de novo binding sites and anchored novel chromatin interactions, and these changes were associated with increased chromatin accessibility and elevated expression of cytotoxic program genes including Tbx21, Ifng, and Klrg1. CTCF was also evicted from its ex-binding sites in naive state, with concomitantly reduced chromatin interactions in effector cells, as observed at memory precursor-associated genes including Il7r, Sell, and Tcf7. Genetic ablation of CTCF indeed diminished cytotoxic gene expression, but paradoxically elevated expression of memory precursor genes. Comparative Hi-C analysis revealed that key memory precursor genes were harbored within insulated neighborhoods demarcated by constitutive CTCF binding, and their induction was likely due to disrupted CTCF-dependent insulation. CTCF thus promotes cytotoxic effector differentiation by integrating local chromatin accessibility control and higher-order genomic reorganization.

The Interplay of Transcription and Genome Topology Programs T Cell Development and Differentiation.

T cells are essential for mounting defense against various pathogens and malignantly transformed cells. Thymic development and peripheral T cell differentiation are highly orchestrated biological processes that require precise gene regulation. Higher-order genome organization on multiple scales, in the form of chromatin loops, topologically associating domains and compartments, provides pivotal control of T cell gene expression. CTCF and the cohesin machinery are ubiquitously expressed architectural proteins responsible for establishing chromatin structures. Recent studies indicate that transcription factors, such as T lineage-defining Tcf1 and TCR-induced Batf, may have intrinsic ability and/or engage CTCF to shape chromatin architecture. In this article, we summarize current knowledge on the dynamic changes in genome topology that underlie normal or leukemic T cell development, CD4+ helper T cell differentiation, and CD8+ cytotoxic T cell functions. The knowledge lays a solid foundation for elucidating the causative link of spatial chromatin configuration to transcriptional and functional output in T cells.

Tcf1-CTCF cooperativity shapes genomic architecture to promote CD8+ T cell homeostasis.

CD8+ T cell homeostasis is maintained by the cytokines IL-7 and IL-15. Here we show that transcription factors Tcf1 and Lef1 were intrinsically required for homeostatic proliferation of CD8+ T cells. Multiomics analyses showed that Tcf1 recruited the genome organizer CTCF and that homeostatic cytokines induced Tcf1-dependent CTCF redistribution in the CD8+ T cell genome. Hi-C coupled with network analyses indicated that Tcf1 and CTCF acted cooperatively to promote chromatin interactions and form highly connected, dynamic interaction hubs in CD8+ T cells before and after cytokine stimulation. Ablating CTCF phenocopied the proliferative defects caused by Tcf1 and Lef1 deficiency. Tcf1 and CTCF controlled a similar set of genes that regulated cell cycle progression and promoted CD8+ T cell homeostatic proliferation in vivo. These findings identified CTCF as a Tcf1 cofactor and uncovered an intricate interplay between Tcf1 and CTCF that modulates the genomic architecture of CD8+ T cells to preserve homeostasis.

Tcf1 preprograms the mobilization of glycolysis in central memory CD8+ T cells during recall responses.

The mechanisms underlying the heightened protection mediated by central memory CD8+ T (TCM) cells remain unclear. Here we show that the transcription factor Tcf1 was required in resting TCM cells to generate secondary effector CD8+ T cells and to clear pathogens during recall responses. Recall stimulation of CD8+ TCM cells caused extensive reprogramming of the transcriptome and chromatin accessibility, leading to rapid induction of glycolytic enzymes, cell cycle regulators and transcriptional regulators, including Id3. This cluster of genes did not require Tcf1 in resting CD8+ TCM cells, but depended on Tcf1 for optimal induction and chromatin opening in recall-stimulated CD8+ TCM cells. Tcf1 bound extensively to these recall-induced gene loci in resting CD8+ TCM cells and mediated chromatin interactions that positioned these genes in architectural proximity with poised enhancers. Thus, Tcf1 preprogramed a transcriptional program that supported the bioenergetic and proliferative needs of CD8+ TCM cells in case of a secondary challenge.

Lysine methyltransferase Kmt2d regulates naive CD8+ T cell activation-induced survival.

Lysine specific methyltransferase 2D (Kmt2d) catalyzes the mono-methylation of histone 3 lysine 4 (H3K4me1) and plays a critical role in regulatory T cell generation via modulating Foxp3 gene expression. Here we report a role of Kmt2d in naïve CD8+ T cell generation and survival. In the absence of Kmt2d, the number of CD8+ T cells, particularly naïve CD8+ T cells (CD62Lhi/CD44lo), in spleen was greatly decreased and in vitro activation-related death significantly increased from Kmt2d fl/flCD4cre+ (KO) compared to Kmt2d fl/flCD4cre- (WT) mice. Furthermore, analyses by ChIPseq, RNAseq, and scRNAseq showed reduced H3K4me1 levels in enhancers and reduced expression of apoptosis-related genes in activated naïve CD8+ T cells in the absence of Kmt2d. Finally, we confirmed the activation-induced death of antigen-specific naïve CD8+ T cells in vivo in Kmt2d KO mice upon challenge with Listeria monocytogenes infection. These findings reveal that Kmt2d regulates activation-induced naïve CD8+ T cell survival via modulating H3K4me1 levels in enhancer regions of apoptosis and immune function-related genes.

Tcf1 and Lef1 provide constant supervision to mature CD8+ T cell identity and function by organizing genomic architecture.

T cell identity is established during thymic development, but how it is maintained in the periphery remains unknown. Here we show that ablating Tcf1 and Lef1 transcription factors in mature CD8+ T cells aberrantly induces genes from non-T cell lineages. Using high-throughput chromosome-conformation-capture sequencing, we demonstrate that Tcf1/Lef1 are important for maintaining three-dimensional genome organization at multiple scales in CD8+ T cells. Comprehensive network analyses coupled with genome-wide profiling of chromatin accessibility and Tcf1 occupancy show the direct impact of Tcf1/Lef1 on the T cell genome is to promote formation of extensively interconnected hubs through enforcing chromatin interaction and accessibility. The integrative mechanisms utilized by Tcf1/Lef1 underlie activation of T cell identity genes and repression of non-T lineage genes, conferring fine control of various T cell functionalities. These findings suggest that Tcf1/Lef1 control global genome organization and help form intricate chromatin-interacting hubs to facilitate promoter-enhancer/silencer contact, hence providing constant supervision of CD8+ T cell identity and function.

Sepsis leads to lasting changes in phenotype and function of memory CD8 T cells.

The global health burden due to sepsis and the associated cytokine storm is substantial. While early intervention has improved survival during the cytokine storm, those that survive can enter a state of chronic immunoparalysis defined by transient lymphopenia and functional deficits of surviving cells. Memory CD8 T cells provide rapid cytolysis and cytokine production following re-encounter with their cognate antigen to promote long-term immunity, and CD8 T cell impairment due to sepsis can pre-dispose individuals to re-infection. While the acute influence of sepsis on memory CD8 T cells has been characterized, if and to what extent pre-existing memory CD8 T cells recover remains unknown. Here, we observed that central memory CD8 T cells (TCM) from septic patients proliferate more than those from healthy individuals. Utilizing LCMV immune mice and a CLP model to induce sepsis, we demonstrated that TCM proliferation is associated with numerical recovery of pathogen-specific memory CD8 T cells following sepsis-induced lymphopenia. This increased proliferation leads to changes in composition of memory CD8 T cell compartment and altered tissue localization. Further, memory CD8 T cells from sepsis survivors have an altered transcriptional profile and chromatin accessibility indicating long-lasting T cell intrinsic changes. The sepsis-induced changes in the composition of the memory CD8 T cell pool and transcriptional landscape culminated in altered T cell function and reduced capacity to control L. monocytogenes infection. Thus, sepsis leads to long-term alterations in memory CD8 T cell phenotype, protective function and localization potentially changing host capacity to respond to re-infection.

A dirty cut, a nasty burn, a severe COVID infection; there are many ways for someone to develop sepsis. This life-threatening condition emerges when the immune system overreacts to a threat and ends up damaging the body. Even when patients survive, they are often left with a partially impaired immune system that cannot adequately protect against microbes and cancer; this is known as immunoparalysis. Memory CD8 T cells, a type of immune cell that is compromised by sepsis, are a long-lived population of cells that 'remember' previous infection or vaccination, and then react faster to prevent the same illness if the person ever encounters the same threat again. Yet it is unclear how exactly sepsis harms the function and representation of memory CD8 T cells, and the immune system in general. Jensen et al. investigated this question, first by showing that sepsis leads to a profound loss of memory CD8 T cells, but that surviving memory CD8 T cells multiply quickly ­ especially a subpopulation known as central memory CD8 T cells ­ to re-establish the memory CD8 T cell population. Since the central memory CD8 T cells proliferate better than the other memory T cells this alters the overall composition of the pool of memory CD8 T cells, with central memory cells becoming overrepresented. Further experiments revealed that this biasing toward central memory T cells, due to sepsis, created long-term changes in the distribution of memory CD8 T cells throughout the body. The way the genetic information of these cells was packaged had also been altered, as well as which genes were switched on or off. Overall, these changes reduced the ability of memory CD8 T cells to control infections. Together, these findings help to understand how immunoparalysis can emerge after sepsis, and what could be done to correct it. These findings could also be applied to other conditions ­ such as COVID-19 ­ which may cause similar long-term changes to the immune system.

Arid1a-Plagl1-Hh signaling is indispensable for differentiation-associated cell cycle arrest of tooth root progenitors.

Chromatin remodelers often show broad expression patterns in multiple cell types yet can elicit cell-specific effects in development and diseases. Arid1a binds DNA and regulates gene expression during tissue development and homeostasis. However, it is unclear how Arid1a achieves its functional specificity in regulating progenitor cells. Using the tooth root as a model, we show that loss of Arid1a impairs the differentiation-associated cell cycle arrest of tooth root progenitors through Hedgehog (Hh) signaling regulation, leading to shortened roots. Our data suggest that Plagl1, as a co-factor, endows Arid1a with its cell-type/spatial functional specificity. Furthermore, we show that loss of Arid1a leads to increased expression of Arid1b, which is also indispensable for odontoblast differentiation but is not involved in regulation of Hh signaling. This study expands our knowledge of the intricate interactions among chromatin remodelers, transcription factors, and signaling molecules during progenitor cell fate determination and lineage commitment.

Interplay of BAF and MLL4 promotes cell type-specific enhancer activation.

Cell type-specific enhancers are activated by coordinated actions of lineage-determining transcription factors (LDTFs) and chromatin regulators. The SWI/SNF chromatin remodeling complex BAF and the histone H3K4 methyltransferase MLL4 (KMT2D) are both implicated in enhancer activation. However, the interplay between BAF and MLL4 in enhancer activation remains unclear. Using adipogenesis as a model system, we identify BAF as the major SWI/SNF complex that colocalizes with MLL4 and LDTFs on active enhancers and is required for cell differentiation. In contrast, the promoter enriched SWI/SNF complex PBAF is dispensable for adipogenesis. By depleting BAF subunits SMARCA4 (BRG1) and SMARCB1 (SNF5) as well as MLL4 in cells, we show that BAF and MLL4 reciprocally regulate each other's binding on active enhancers before and during adipogenesis. By focusing on enhancer activation by the adipogenic pioneer transcription factor C/EBPß without inducing cell differentiation, we provide direct evidence for an interdependent relationship between BAF and MLL4 in activating cell type-specific enhancers. Together, these findings reveal a positive feedback between BAF and MLL4 in promoting LDTF-dependent activation of cell type-specific enhancers.

Evidence of Superior and Inferior Sinoatrial Nodes in the Mammalian Heart.

OBJECTIVES: This study sought to investigate the shift of leading pacemaker locations in healthy and failing mammalian hearts over the entire range of physiological heart rates (HRs), and to molecularly characterize spatial regions of spontaneous activity. BACKGROUND: A normal heartbeat originates as an action potential in a group of pacemaker cells known as the sinoatrial node (SAN), located near the superior vena cava. HRs and the anatomical site of origin of pacemaker activity in the adult heart are known to dynamically change in response to various physiological inputs, yet the mechanism of this pacemaker shift is not well understood. METHODS: Optical mapping was applied to ex vivo rat and human isolated right atrial tissues, and HRs were modulated with acetylcholine and isoproterenol. RNA sequencing was performed on tissue areas that elicited spontaneous activity, and comparisons were made to neighboring myocardial tissues. RESULTS: Functional and molecular evidence identified and confirmed the presence of 2 competing right atrial pacemakers localized near the superior vena cava and the inferior vena cava-the superior SAN (sSAN) and inferior SAN (iSAN), respectively-which preferentially control the fast and slow HRs. Both of these regions were evident in non-failing rat and human hearts and maintained spontaneous activity in the rat heart when physically separated from one another. Molecular analysis of these 2 pacemaker regions revealed unique but similar transcriptional profiles, suggesting iSAN dominance when the sSAN is silent. CONCLUSIONS: The presence of 2 spatially distinct dominant pacemakers, sSAN and iSAN, in the mammalian heart clarifies previous identification of migrating pacemakers and corresponding changes in P-wave morphology in mammalian species.

Transcriptome-wide stability analysis uncovers LARP4-mediated NFκB1 mRNA stabilization during T cell activation.

T cell activation is a well-established model for studying cellular responses to exogenous stimulation. Motivated by our previous finding that intron retention (IR) could lead to transcript instability, in this study, we performed BruChase-Seq to experimentally monitor the expression dynamics of nascent transcripts in resting and activated CD4+ T cells. Computational modeling was then applied to quantify the stability of spliced and intron-retained transcripts on a genome-wide scale. Beyond substantiating that intron-retained transcripts were considerably less stable than spliced transcripts, we found a global stabilization of spliced mRNAs upon T cell activation, although the stability of intron-retained transcripts remained relatively constant. In addition, we identified that La-related protein 4 (LARP4), an RNA-binding protein (RBP) known to enhance mRNA stability, was involved in T cell activation-dependent mRNA stabilization. Knocking out Larp4 in mice destabilized Nfκb1 mRNAs and reduced secretion of interleukin-2 (IL2) and interferon-gamma (IFNγ), two factors critical for T cell proliferation and function. We propose that coordination between splicing regulation and mRNA stability may provide a novel paradigm to control spatiotemporal gene expression during T cell activation.

HDAC10 Regulates Cancer Stem-Like Cell Properties in KRAS-Driven Lung Adenocarcinoma.

Activation of oncogenic KRAS is the most common driving event in lung adenocarcinoma development. Despite the existing rationale for targeting activated KRAS and its downstream effectors, the failure of clinical trials to date indicates that the mechanism of KRAS-driven malignancy remains poorly understood. Here we report that histone deacetylase 10 (HDAC10) might function as a putative tumor suppressor in mice carrying a spontaneously activated oncogenic Kras allele. Hdac10 deletion accelerated KRAS-driven early-onset lung adenocarcinomas, increased macrophage infiltration in the tumor microenvironment, and shortened survival time in mice. Highly tumorigenic and stem-like lung adenocarcinoma cells were increased in Hdac10-deleted tumors compared with Hdac10 wild-type tumors. HDAC10 regulated the stem-like properties of KRAS-expressing tumor cells by targeting SOX9. Expression of SOX9 was significantly increased in Hdac10-deleted tumor cells and depletion of SOX9 in Hdac10 knockout (KO) lung adenocarcinoma cells inhibited growth of tumorspheres. The genes associated with TGFß pathway were enriched in Hdac10 KO tumor cells, and activation of TGFß signaling contributed to SOX9 induction in Hdac10 KO lung adenocarcinoma cells. Overall, our study evaluates the functions and mechanisms of action of HDAC10 in lung carcinogenesis that will inform the rationale for targeting its related regulatory signaling as an anticancer strategy. SIGNIFICANCE: These findings linking HDAC10 and lung tumorigenesis identify potential novel strategies for targeting HDAC10 as a treatment for lung cancer.

RECOGNICER: A coarse-graining approach for identifying broad domains from ChIP-seq data.

BACKGROUND: Histone modifications are major factors that define chromatin states and have functions in regulating gene expression in eukaryotic cells. Chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq) technique has been widely used for profiling the genome-wide distribution of chromatin-associating protein factors. Some histone modifications, such as H3K27me3 and H3K9me3, usually mark broad domains in the genome ranging from kilobases (kb) to megabases (Mb) long, resulting in diffuse patterns in the ChIP-seq data that are challenging for signal separation. While most existing ChIP-seq peak-calling algorithms are based on local statistical models without account of multi-scale features, a principled method to identify scale-free board domains has been lacking. METHODS: Here we present RECOGNICER (Recursive coarse-graining identification for ChIP-seq enriched regions), a computational method for identifying ChIP-seq enriched domains on a large range of scales. The algorithm is based on a coarse-graining approach, which uses recursive block transformations to determine spatial clustering of local enriched elements across multiple length scales. RESULTS: We apply RECOGNICER to call H3K27me3 domains from ChIP-seq data, and validate the results based on H3K27me3's association with repressive gene expression. We show that RECOGNICER outperforms existing ChIP-seq broad domain calling tools in identifying more whole domains than separated pieces. CONCLUSION: RECOGNICER can be a useful bioinformatics tool for next-generation sequencing data analysis in epigenomics research.